Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 25 February 2008 



(MN WF$i style file v2.2) 



The Impact of Stellar Oscillations on Doppler Velocity 
Planet Searches 



oo 
O 

o 

(N 

X> 
<D 
IX, 

m 

(N 

Oh 
6 

o3 



(N 
> 
OO 

m 

m 

o 
o 



X 



S. J. O'Toole, 1 C. G. Tinney, 1 * 2 H. R. A. Jones 3 

1 Anglo -Australian Observatory, PO Box 296, Epping 1710, Australia 

2 Department of Astrophysics, School of Physics, University of NSW, 2052, Australia 

3 Centre for Astrophysical Research, University of Hertfordshire, Hatfield, AL10 9AB, UK 



25 February 2008 



ABSTRACT 

We present a quantitative investigation of the effect of stellar oscillations on Doppler 
velocity planet searches. Using data from four asteroseismological observation cam- 
paigns, we find a power law relationship between the noise impact of these oscillations 
on Doppler velocities and both the luminosity-to-mass of the target stars, and observed 
integration times. Including the impact of oscillation jitter should improve the quality 
of Keplerian fits to Doppler velocity data. The scale of the effect these oscillations have 
on Doppler velocity measurements is smaller than that produced by stellar activity, 
but is most significant for giant and subgiant stars, and at short integration times (i.e. 
less than a few minutes). Such short observation times tend to be used only for very 
bright stars. However, since it is these very same stars that tend to be targeted for 
the highest precision observations, as planet searches probe to lower and lower planet 
masses, oscillation noise for these stars can be significant and needs to be accounted 
for in observing strategies. 
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1 INTRODUCTION 



Uncertainties in high-precision Doppler velocity measure- 
ments are now reaching ~lms"' on a regular basis (e.g. 
lO'Toole etaLll2007l ; IPepe et al.ll2007f ). The science drivers 
behind the quest for ever-greater precision are the detection 
of Earth-mass planets in short-period orbits and Solar Sys- 
tem analogues, as well as very low amplitude stellar oscilla- 
tions. Understanding both increasingly subtle instrumental 
variations, and intrinsic stellar variability, is now more im- 
portant than ever. 

The term jitter has been coined to describe the 
noise imposed on precision radial velocity programs by 
a star's intrinsic instab ility and has b e en in vestigated by 
ISaar fc Donahue! (|l997h and lSaar etaD l| 19981 ). with the ef- 
fect particularly detrimental on giant stars. Until now, only 
the stellar acti vity componen t of jitter has been investigated 
quantitatively. IWrightl (|2005h examined the jitter of targets 
in the Lick and Keck Planet Searches and derived an em- 
pirical relationship based on colour, activity and absolute 
magnitude. This has allowed the inclusion of an additional 
term to the measurement uncertainties used when fitting 
Keplerians to velocity data, and therefore better modelling 
of the scatter about these fits. 

The dominant oscillations in the solar-like stars are non- 
radial p-modes stochastically excited by turbulent convec- 



tion. At least 12 stars have been obser ved by various groups 
looki ng for solar-like oscillations (cf. Bed ding fc Kieldsenl 
2006). The amplitudes of these velocity variations depend 
on the stellar luminosity and mass and can be scaled from 
the Sun, at least for objects that have not evolv ed too far off 
the main sequence (jKieldsen fc Bedding! 1 1995T ). These suc- 
cessful detections of solar-like oscillations have been in no 
small part due to advances in precision radial velocity tech- 
niques made by planet search teams; however, little quan- 
titative work has been done to examine the impact of the 
osci llations on the detec tion of extra-solar planets. 

iTinnev et all (|2005l ) discussed asteroseismology as noise 
in the context of Doppler velocity planet searches. They sug- 
gested that exposure times of integer multiples of the peak 
oscillation period s could lower j itter by 1-2 ms . This was 
also discussed bv lMavor et al.l (|2003l ). who argued that ex- 
posure times of around 15 minutes would average out oscilla- 
tions. To date these suggestions have not been quantitatively 
investigated. 



2 OBSERVATIONS 

To quantify the impact of p-mode oscillations on Doppler 
planet searches as a function of both observing strategy and 
stellar properties we have analysed data from four of the 



2 S. J. 'Toole et al 




3. 10.0 : 



* X □ 



X « Ceil R 
□ a Cen A 
A (3 Hyi 
o f Ind 



1709.0 1709.1 1709. 'i 

Time (days) 



1709.3 1709.4 



Window size (mi] 



Figure 1. Moving average with 600 s window for a set of obser- 
vations of a Cen A and /3 Hyi. 



Figure 2. The total velocity range in an observation data set for 
each star, as a function of the moving average window size. 



published programs of high time-resolution observations ob- 
tained with the University College London Echelle Spectro- 
graph (UCLES) that detecte d and analysed so lar-like os- 
cillatio ns in a Cen A (G2V ; iButler et all |2004). a Cen B 
(K1V ; iKieldsen et al.|[2005l), j3 Hyi (G2IV; iBedding et all 
l200ll ) and v Indi fGOIVi lBedding et al.ll2006l ). We note that 
while several of these papers used data obtained with the 
CORALIE and UVES spectrographs, we have access only 
to the UCLES data, so these other data are not analysed 
here. 

The observations examined in this paper are described 
in detail in the references giv en above and are a lmost the 
same as used for the AAPS IButler et all 1120011 ). Briefly, 
the data were taken using UCLES mounted at the coude 
focus of the Anglo- Australian Telescope ( AAT) . An iodine 
absorption cell is placed in the beam, imprinting a forest of 
molecular iodine absorption lines onto the stellar spectrum. 
These lines are used as a wavelengt h reference to derive high- 
precision velocities as described in IButler et al.l |T996). The 
integration times for each object depend on a number of 
factors including its brightness, its expected dominant os- 
cillation period, P ma x and current weather conditions. The 
data analysed in this paper is exactly the same as that used 
for the asteroseismological analyses described in the refer- 
ences above, including the removal of long-term drifts in the 
velocity time series for all stars except v Indi (where such a 
correction was not found necessary). 



3 THE EFFECT OF OSCILLATIONS 

Observed stellar radial velocity curves show variations from 
several sources: those intrinsic to the star; reflex motion due 
to a companion; and changes and drifts in the instrumental 
setup. This paper investigates the observational effect of the 
first of these - in particular, stellar oscillations - on the 
second. 



3.1 Oscillations as noise 

To look at the impact of exposure time on the velocity varia- 
tions that p-mode oscillations produce in these stars, we have 



calculated a moving average of our high time-resolution as- 
teroseismology data sets. The window of the moving average 
was set to typical AAPS exposure times (5, 10, 15, 20, 30, 
45, and 60 minutes) as well as several longer times to exam- 
ine the effects of averaging an entire night's asteroseismology 
observations (90, 120, 180, 300, 450 and 600 minutes). An 
example of averaged data over time for a Cen A and /3 Hyi is 
shown in Figure[T]with a window of 10 minutes. The times- 
tamp is set to the midtime of the observations in the window. 
Ten minute exposures sample almost a cycle and a half of 
the dominant periodicity in a Cen A, but significantly less 
than a cycle for f3 Hyi (taking the dominant period as the 
individual mode with the highest amplitude). In both cases, 
there is still significant scatter in the time series when the 
total integration time spans ten minutes. 

If we measure the total velocity range spanned by the 
velocity extrema for each star and each windowing time 
listed above, we produce the results shown in Figure[2] These 
extrema are measured from all nights in the time series. If 
we consider the 10 minute window for a Cen A shown in Fig- 
ure [TJ we see that the star varies by up to around 5ms _1 . 
This number decreases conside rably the more cycle s we "in- 
tegrate" over, as suggested by I Thine v et alj (2005). 

The size of the variations should vary with spec- 
tral type, since oscillation amplitudes are dependent on 
stellar luminosity and m ass. Based on equation 7 of 
iKield sen fc Beddind (| 19951 ). earlier- type stars should be 
more affected than later-type stars, and subgiants can be 
expected to show the largest effects. From Figure [2] we see 
that this is the case and that a Cen A has higher scatter due 
to oscillations than the later-type star a Cen B - particu- 
larly at short total integration times. Subgiants show larger 
scatter due to oscillations than main-sequence stars, with 
the more evolved subgiant v Ind showing more scatter than 
the less evolved /3 Hyi. We note that while the scatter is 
problematic for planet searches, the oscillations which cause 
it allow precise determination of the stellar mass, in turn 
leading to more precise planetary masses. 
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Figure 3. (a) Ninety-five percent confidence ranges (n osc ) for each star as a function of simulated integration time (tint), (b) Data from 
the panel (a) normalised to produce n' OBC and averaged. Uncertainties represent the standard deviation; the standard error of the mean 
is a factor of two smaller. 



3.2 Confidence limits and an empirical relation 

While the velocity range data of Figure [2] reveal the im- 
portance of stellar oscillations to precision Doppler planet 
searches, what is really required is an understanding of their 
statistical impact. That is, how can they be modelled as a 
source of Doppler noise for planet searches? 

Doppler programs typically model their noise sources as 
Gaussian distributions, and in an ideal world, the "noise im- 
pact" of these stellar oscillations would be parameterised in 
the same way. However, as even a cursory gla nce at a stel- 
lar o scillation power spectrum indicates (e.g. iButler et al.l 
120041 ). they are typically not a source of Gaussian noise. 
Here we are not referring to the underlying noise of the ob- 
servations in the absence of oscillations, rather to the con- 
tribution of the oscillations themselves. We have therefore 
parameterised the impact of oscillations using "95% confi- 
dence velocities" - i.e. the velocity range, n osc , within which 
95% of the measured velocities for a given target would lie 
for a given simulated integration time, tint El 

Figure[2]suggests there exists a quantifiable relationship 
between the jitter due to solar-like oscillations and total in- 
tegration time. To examine this Figure [3j a) shows ^osc as a 
function of tint for all four stars. There appears to be a con- 
sistent trend in each case, especially at periods of 5-60 min. 
This is not all that surprising, since while the detailed form 
of the stellar oscillation power spectra in these stars (e . g. fine 
structure splitting between modes; iKieldsen et all 120051 ) 
will depend on the details of their interior structure, the 
overall envelope of their power spectra (which is what we 



1 Note that if the asteroseismological power spectrum were Gaus- 
sian, a 95% confidence velocity corresponds to a velocity range of 
±1.96 cr about the mean velocity, or that tJ OBC ~ n osc /4. 



sample in these observations) are very similar. We have nor- 
malised each star's n osc values at ti nt = 30 min (which is a 
typical longest exposure time in the AAPS for our very high- 
est precision targets) to produce n' osc , which we then average 
over all four stars and plot in Figure |3jb); the uncertainties 
are a simple standard deviation. We have modelled this as 
two linear trends with a break point at 35 minutes. We find 



3.20- 1.51 log 10 



for tint < 35 min 



2.37 - 0.99 log 10 tint for 35 < t in t < 180 min 



(1) 



It is clear from Figure [3] that this relationship breaks down 
above tint ~ 180 min. This is not surprising, as almost all of 
the power above this time-scale has been extracted by the 
high-pass filtering of the asteroseismology data to remove 
long-term drifts. 

We expect "oscillation noise" to depend on the ratio of 
stellar luminosity t o mass in a similar manne r to the oscil- 
lations themselves |Kield sen fc Beddindll995T) . To examine 
this, we plot n osc /4 in FigureUfor each star as a function of 
its luminosity-mass ratio (L/M), at a range of simulated ex- 
posure times. (n osc /4 is plotted as this form is more directly 
comparable with 1-cr noise estimates from other sources.) 

That n osc /4 has approximately the same slope for each 
exposure time on a logarithmic scale suggests a power law 
relationship exists with L/M. We fit a power law at ti nt = 
30 min and find 



n' OEC x 10 311 (L/M)° 



for tint < 180 min 



(2) 



where n' osc is given by equation [T] The fit is overplotted in 
Figure [4] along with power laws scaled to the other values of 
tint shown using the appropriate n' osc . 

Also shown in Figure [3] are the residual root- 
mean-square (RMS ) values for each known planet from 
IButler et ail ((2006). The RMS represents the average de- 
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viation from a perfect fit to the radial velocity data and is 
made up of several components, including the various forms 
of jitter, instrumental variations and data quality. It is also 
worth noting here that low-mass undetected planets are also 
a source of noise. 

The impact of oscillation jitter can be seen (in general) 
to skirt the lower edge of the observed exoplanetary RMS 
values, though for very short exposure times, or evolved 
stars, the noise impact reaches amplitudes of several ms -1 , 
where it clearly becomes significant and of concern. 
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4 DISCUSSION 

We have demonstrated that the impact of p-mode oscilla- 
tions on low-amplitude planet searches can be quantified. 
Unlike the stellar activity jitter however, oscillation jitter 
is dependent on the length of time spent observing a tar- 
get at any given epoch, so it therefore affects the observing 
strategies for planet hunting. 



Figure 4. Oscillation jitter (n OS c/4 PS ct sc) of the UCLES as- 
teroseismology targets as a function of log 10 (L/M) for various 
simulated integration times: 5 minutes (crosses); 10 min (trian- 
gles); 30 min (circles); 60 min (squares); 120 min (asterisks); 300 
min (plus signs). Overplotted are the power laws dervied from 
Equation [2] Final l y, the residual RMS values for known planets 
from I Butler et all ||2006| . small diamonds) are also shown. 



4.1 Implications for Observing Strategies 

Observing strategies should be tailored to minimise the im- 
pact of oscillations. The key factors to consider are the tar- 
get's L/M value and the resulting total integration time 
needed to lower 7i osc to a point where it drops below the 
desired photon-counting signal-to-noise ratio (SNR) require- 
ment. It is very unlikely that lengthy asteroseismology cam- 
paigns will be staged for even a small fraction of the stars; 
however, one can use Equations [T] and [2] to optimise inte- 
gration times. More evolved stars will have a higher L/M 
and will be more affected by oscillation jitter, therefore re- 
quiring longer integration times, re gardless of the brightness 
of the object. iHekker et a l. (2006) discussed radial velocity 
variations in giant stars and suggest that variations could 
be even larger than predicted here for these objects. 

Consider, for example, the bright, slightly evolved star fj, 
Ara (HD 160691; V=5.12; L = 1.75L Q ; M = 1.15M Q ) which 
has been the subject of numerous planet discovery papers, 
and has been claimed to host up to 4 planets, the smallest 
producing velocity amplitudes as low as 3ms -1 . This star 
is very bright, and so requires integration times of only a 
few minutes to reach a SNR sufficient to achieve ~ 1ms" 1 
precision or better with the AAT and UCLES. At t int = 1 & 
5min, we find n osc /4 = 0.58, 0.39ms -1 - sufficient to make 
constraining the innermost planet in this system difficult or 
impossible. Observations of ij n t > 15 min are required to 
reduce n OBC /4 below 0.26 ms -1 , and of > 60 min to reduce 
n osc /4 below 0.11ms -1 . 

As seen above, the scale of the effect these oscillations 
have on Doppler velocity measurements is smaller than that 
produced by stellar activity, but is most significant for giant 
and sub-giant stars, and at short integration times (i.e. less 
than a few minutes). Such short observations times tend 
to be used only for very bright stars. However, as planet 
searches target lower and lower masses, it is these very same 
stars that tend to be targeted for the highest precision obser- 
vations. So oscillation noise for these stars can be important 
and needs to be accounted for in observing strategies. 



4.2 Avenues of Further Investigation 

Apart from solar-like oscillations and stellar activity, what 
other stellar noise sources remain to be quantified? Two 
sources of potential noise are stellar a ctivity cycles a nd stel- 
lar convection. The activity metric of I Wright] l|2005h is use- 
ful and has been widely adopted; however, it does not in- 
clud e a time-varying com ponent, which is certainly present 
fe.g. lMetcalfe et al.| [2007l. The timescales of these variations 
though are much longer than oscillations - the order of years 
rather than minutes. Incorporating a time dependence into 
jitter measurements is worthy of investigation, especially 
since stellar activity timescales are similar to Jupiter's or- 
bital period. 

Sun-like stars have large convective cells or granules 
where material is dredged up from lower layers and mixed to 
the surface. The process involves many random surface mo- 
tions that even when averaged over time are almost certainly 
large enough to affect planet search measu rements. Vari - 
ations are expected to be around 1-2 ms" _1 <lpravinslll999r) 
and m ay have already been observed in p Ara ( Bouchv et al.l 
120051 ). Granule lifetimes are typically tens of minutes which 
is similar to the exposure times used in planet searches. De- 
spite these characteristics, the impact of convection has not 
been quantitatively investigated. 

These two effects will form the focus of our ongoing 
investigation of the noise sources limiting precision Doppler 
planet searches. 



5 CONCLUSIONS 

We have shown that the noise impact of stellar oscillations 
on precision Doppler velocities obtained in the search for 
extra-solar planets can be significant in some circumstances. 
We have used asteroseismological data sets to derive rela- 
tions which quantify that impact as a function of integra- 
tion time, and stellar luminosity-to-mass ratio. These can 
be used to improve the quality of Keplerian fits to planet 
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search data. But most importantly, these relations can drive 
observing strategies in the search for the lowest-mass planets 
around bright and evolved stars. 
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